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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  MEMORANDUM  1406 


ON  THE  USE  OF  THE  HARMONIC  LINEARIZATION  METHOD  IN  THE 

AUTOMATIC  CONTROL  THEORY* 

By  E.  P.  Popov 

The  method  of  harmonic  linearization  (harmonic  balance),  first 
proposed  by  N.  M.  Krylov  and  N.  N.  Bogolyubov  (ref .  l)  for  the  approxi- 
mate investigation  of  nonlinear  vibrations,  has  been  developed  and  re- 
eived  wide  practical  application  to  problems  in  the  theory  of  automatic 
control  (refs.  3  to  6) .   Recently,  some  doubt  has  been  expressed  on  the 
legitimacy  of  application  of  the  method  to  these  problems,  and  assertions 
were  made  on  the  absence  in  them  of  a  small  parameter  of  any  kind.   Never- 
theless, the  method  gives  practical,  acceptable  results  and  is  a  simple 
and  powerful  means  in  engineering  computatations .   Hence,  the  importance 
of  questions  arises  as  to  its  justification.   The  underlying  principle  of 
the  method  is  the  replacement  of  the  given  nonlinear  equation  by  a  linear 
equation.   In  establishing  the  method,  a  small  parameter  is  considered 
whose  presence  makes  it  possible  to  speak,  with  some  degree  of  approxi- 
mation, of  the  solution  of  this  new  equation  to  the  solution  of  the  given 
nonlinear  equation.   In  an  article  by  the  author  (ref.  7),  certain  con- 
siderations were  given  on  the  presence  of  the  small  parameter,  but  this 
question  has  not  as  yet  received  a  final  answer.   In  the  present  report, 
a  somewhat  different  approach  to  the  problem  is  applied  that  permits: 
(a)  establishing,  in  the  clearest  manner,  the  form  of  the  presence  of 
the  small  parameter  in  nonlinear  problems  of  control  theory,  solvable  by 
the  method  of  harmonic  linearization;  (b)  connecting  it  with  previous 
intuitive  physical  concepts  (with  the  "filter  property")  and  extending 
the  class  of  problems  possessing  this  property;  and  (c)  discussing  various 
generalizations  of  the  method. 

The  free  motion  (transition  process  and  autovibration)  for  a  very 
wide  class  of  nonlinear  systems  of  automatic  control  (ref.  7)  are 
described  by  differential  equations  of  the  form 


Q(p)x  +  R(p)F(x,px)  =0   (p  =  ^)  (1) 


* "K  voprosu  o  primenenii  metoda  harmonicheskoi  linearizatsii  v  teorii 
regulirovaniya."  Doklady  Akademii  Nauk  (SSSR),  vol.  106,  no.  2,  1956,  pp. 
211-214. 
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where   Q(p)  and  R(p)  are  polynomials  of  any  degree,  of  which  some  re- 
quired properties  will  be  established  in  the  following  paragraphs,  and 
F(x,px)  is  a  given  nonlinear  function  possible  only  with  respect  to  as- 
sumptions of  the  most  general  character.   However,  in  problems  of  the 
theory  of  control,  no  assumptions  must  be  made  as  to  the  smallness  of 
the  nonlinear  function  F(x,px)  or  to  its  small  difference  from  a  linear 
function.   In  order  to  render  explicit  the  form  in  which  it  would  be 
possible  to  write  a  small  parameter  in  equation  (l),  we  shall  proceed  as 
follows . 

Let  equation  (l)  have  a  periodic  solution  or  a  solution  approximately 
periodic  differing  slightly  from  the  sinusoidal.  We  write  this  solution 
in  the  form 

X  =  X*  +  ey(t)    X  =  a  sin  oot  (2) 

where  e   denotes  a  small  parameter,  and  y(t)  denotes  an  unknown  bounded 
function  of  time.   In  the  case  of  the  existence  of  a  periodic  solution, 
we  write 

00 

ey/t)    =  e^  a^  sin(koDt  +  <!>^)  (3) 

k=2 

We  represent  the  given  nonlinear  function  F(x,px)  in  the  form 

F(x,p)  =  F(x*,px*)  +  [F(x*  +  £y,  px*  +  epy)  -F(x*,px*)]    (4) 

Expanding  the  two  components  separately  in  a  Fourier  series,  we  obtain 


F(x,px)  =  q^  +  (A  +  ^  p)  sin  oot  +  2_.  F,  +  &Y2 
\      /         k=2       k=0 


\  (5) 


k 


where  q^.  A,  and  B  are  the  coefficients  of  the  initial  term,  the  sine, 

and  cosine  terms,  rexpectively ,  the  cosine  being  replaced  by  —  p  sin  oot 

of  the  expansion  of  the  function  F(x*,px*)  in  a  Fourier  series.  ZF-^ 
denotes  all  the  higher  harmonics  of  the  expansion  of  F(x*,px*)  in  a 
Fourier  series  (they  must  not  be  considered  small,  since  the  nonlinearity 
is  not  small)  where  we  write 

F^  =   bj^  sin  (koit  +^^)  (k  =  2,  3,  .  .  .,  -)  (6) 

e2§,   denotes  all  the  terms  of  the  expansion  in  a  Fourier  series  of  all 
the  expressions  shown  in  brackets  in  formula  (4) .   This  entire  expression 
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is  written  with  a  small  parameter  which,  according  to  equation  (4),  is 
small  if  the  derivatives  Sf/Sx  and  SF/Spx  are  finite.   This  expression 
is  also  computed  as  small  in  the  case  of  certain  discontinuous  nonlinear 
characteristics  (e.g.,  the  Raleigh  type  where  the  preceding  derivatives 
at  the  points  of  discontinuity  are  delta  functions) .  We  may  write 


e* 


^  =  ecj^  sin(tojat  +  ^^)  (k  =  0,  1,  2, 


,-) 


(7) 


We  substitute  equations  (2)  and  equation  (4)  in  the  given  equation  (l) , 
so  that 


Q(p)x*  +  Q(p)ey  + 

I  \  OP  OP 

R(p)qO  +  R(P)U  +  ^  p)  sin  oDt  +  R(p)e^  F^  +   R(p)£^<i>j^  =  0 


(8) 


k=2 


k=0 


Since  the  equation  must  be  satisfied  identically,  we  separately 
equate  to  zero  all  coefficients  with  the  same  order  harmonic.   We  note 
that  formula  (8) ,  in  the  case  of  the  existence  of  a  periodic  solution, 
is  exact . 

From  the  equating  of  the  zeroth  harmonics  of  equation  (8),  there  is 
obtained  with  an  accuracy  up  to  e   the  relation 


(9) 


which  is  a  certain  general  requirement  for  F(x,px) . 

From  the  equating  of  the  first  harmoics  of  equation  (8) ,  taking 
account  of  equations  (2)  and  equation  (7),  we  have 


a  sin  oot  = 


-   A^ 


A'^  +  B^ 


R(iao) 


Q(icjo) 


sin  (out  +  y  +  p)  - 


ECn 


R(icu) 


Q(ia)) 


sm 


(oot 


+  ^3_  + 


P) 


(10) 


where  x     ^^^     P  s-^e  arguments  of  the  expressions  A  +  iB  and 
R(icD)/Q(iaL))  .   On  the  basis  of  the  exact  equation  (lO),  we  obtain  the 
following  approximation  (with  an  accuracy  up  to  e): 


-/- 


A^  +  B^ 


R(iGD) 


Q(ico) 


r  +  P 


(11) 
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It  is  here  assumed  that  the  polynomial  Q(p)  in  equation  (l)  does  not 
have  purely  imaginary  roots . 

From  the  equating  of  the  higher  harmonics  of  equation  (8),  consider- 
ing equations  (3),  (s),  and  (7),  we  have 


■.a^  sin(ka3t   +  'Pv) 


R(ika), 


Q(ikjD 


sin(ka)t  +  \  +  P^^)-ecj^ 


Rfikoj) 


Q(ika)) 


sin(koDt   +  *i^  +  Pi^)      (12) 


where  ^-^     denotes  the  argument  of  the  expression  R(ikJo)/Q(il2Jo)  . 

It  is  thus  seen  that  if  bv_  is  not  small,  the  magnitude 
I  R(ika))/Q(ikco)|   should  be  of  the  order  of  e.   The  last  component  in 
equation  (12)  will  then  be  of  the  order  e.^ .      From  the  exact  equation 
(12) ,  we  obtain  the  following  approximate  equation  (with  an  accuracy  up 
to  e)  : 


e  a^  =  b^ 


R(ika)) 


Q(ik£ja) 


\  +  Pk  '  '^ 


(13) 


Comparing  formulas  (l3)  and  (ll)  we  see  that,  for  example,  the  wish 
to  have  in  the  solution  (see  eqs.  (2)  and  (3)) 


(eaj'« 


(14) 


leads  to  the  need  of  satisfying,  in  the  given  equation  (l),  the  following 
requirement: 


oo 

L 

k=2 


R(ikcD) 
Q(ikGD) 


■<  (a2  +  b2) 


R(ia)) 


Q(icD) 


(15) 


where  its  satisfying  in  the  concrete  system  can  be  checked  after  cd  is 
obtained.   The  degree  of  Q(p)  should,  in  any  case,  be  higher  than  that 
of  R(p) .   A  particular  case  of  the  general  expression  (l5)  is  intuitively 
the  earlier  introduced  "filter  property"  of  the  linear  part. 

Thus,  condition  (l5)  has  been  obtained  and  must  be  satisfied  by  the 
coefficients  of  the  given  differential  equation  (l)  in  order  that  a  pe- 
riodic solution,  if  it  exists,  may  be  approximately  determined  in  the 
form  of  sinusoids  in  the  presence  of  a  "strong"  nonlinearity  of  F(x,px) . 
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The  equation  for  its  approximate  determination  according  to  equation  (s) 
with  the  substitution  of  sin  oot  =  x'Va,  assumes  the  form 

[q(p)  +  R(p)  (q  +  1^  p)]  X*  =  0  (16) 


where 


'2jt 
^  ~  a  ~  T^  I   -^(^  ^^^  ^>    ^^   ^°^  u)sin  u  du 


F(a  sin  u,  aoo  cos  u)cos  u  du  (l7) 

The  replacement  of  equation  (l)  by  equation  (16)  with  its  subsequent 
investigation  by  the  linear  methods  is  called  harmonic  linearization. 
This  is  formally  equivalent  to  following  the  mode  of  writing  the  initial 
equation  (l)  with  a  small  parameter,  so  that 

Q(p)x  +  R(p) fq  +  _i  d  X  +  ef (x,p)  =  0  (18) 

where,  according  to  equation  (8),  we  may  write 

ef(x,p)  =  R(p) 


22 ^k  +  ^ S\  - (<i  +  of  p) ^y 

k=2      k=0    ^      ^      J 


where  the  first  terra  of  the  three  terms  shown  in  brackets,  taken  sepa- 
rately, is  not  small.   The  smallness  of  the  function  ef(x,p)  is  acquired 
only  with  the  factor  R(p)  owing  to  the  property  of  equation  (15) . 

There  has  thus  been  found  the  form  of  the  presence  of  the  small 
parameter   e  in  the  equations  of  nonlinear  systems  of  automatic  control 
required  for  the  application  of  the  method  of  harmonic  linearization 
(this  also  refers  to  the  first  approximation  of  the  method  of  small  pa- 
rameter which,  according  to  reference  7,  coincides  with  the  given  method) . 
The  writing  of  equation  (l)  with  small  parameter  in  the  form  of  equation 
(18)  is  valid  in  the  region  of  the  existence  of  equation  (2). 

For  definiteness  we  shall  assume  that  the  polynomials  R(p)  and 
Q(p)  are  such  that  the  characteristic  equation  (eq.  (16))  has  all  positive 
coefficients,  and  that  the  quotient  of  the  division  of  the  entire  left 
side  of  equation  (16)  by  p^  +  oo^  satisfies  the  criterion  of  Hurwitz. 
It  may  then  be  said  that  the  given  nonlinear  system  (eq.  (l) )  for  the 
existence  of  the  periodic  solution  (eq.  (2))  is  close  to  the  equivalent 
linear  system  which  is  on  the  boundary  of  stability,  except  for  the 
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circumstance  that  the  linear  system  is  different  for  various  periodic 
solutions  (because  of  a  change  in  the  coefficients  q  and  q,   from  one 
periodic  solution  to  another^  since  q  and  q-,   depend  on  the  amplitude 
a  and  the  frequency  oo  of  the  required  solution).   Thus,  evidently,  the 
equation  of  the  first  approximation  (eq.  (l6))  is  nonlinear  if  one  speaks 
of  the  combination  of  all  possible  periodic  solutions  for  different  values 
of  the  coefficients  of  the  polynomials  R  and  Q  (i.e.,  for  different 
parameters  of  the  system) ,  and  is  linear  only  for  each  given  periodic 
solution  (along  the  given  solution) .   Such  is  the  property  of  the  method 
of  harmonic  linearization  in  its  application  to  systems  with  "strong" 
nonlinearity . 

By  writing  the  given  equation  (l)  in  the  form  of  equation  (8)  all 
the  variants  of  the  method  of  harmonic  balance  in  control  theory  and  all 
generalizations  of  the  method  also  become  easily  understandable. 

A  similar  reasoning  may  also  be  developed  for  the  application  of  the 
method  of  N.  N.  Bogolyubov  (ref .  2)  to  the  investigation  of  transient 
vibrational  processes  in  certain  nonlinear  systems  of  control. 
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